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i Abstract 

We present one-loop corrections to the Higgs boson masses in the MSSM with Non- 
Minimal Flavor Violation. The flavor violation is generated from the hypothesis of 
general flavor mixing in the squark mass matrices, and these are parametrized by a 
complete set of 5^ Y (X, Y = L,R; i, j = t,c,u or b, s, d). We calculate the corrections 
to the Higgs masses in terms of these 5^ Y taking into account all relevant restrictions 
from B-physics data. This includes constraints from BR(i3 — > X s j), BR(i? s — > 
and AMb s - After taking into account these constraints we find sizable corrections to 
the Higgs boson masses, in the case of the lightest MSSM Higgs boson mass exceeding 
tens of GeV. These corrections are found mainly for the low tan (3 case. In the case 
of a Higgs boson mass measurement these corrections might be used to set further 
constraints on d£ Y ■ 
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1 Introduction 



The Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) [T] with two 
scalar doublets accommodates five physical Higgs bosons. In lowest order these are the light 
and heavy CP-even h and H, the CP-odd A, and the charged Higgs bosons H ± . The Higgs 
sector of the MSSM can be expressed at lowest order in terms of the gauge couplings, Ma and 
tan/3 = vz/vi, the ratio of the two vacuum expectation values. All other masses and mixing 
angles can therefore be predicted. Higher-order contributions can give large corrections to 
the tree-level relations (see e.g. Ref. [2] for reviews). 

The MSSM predicts scalar partners for each fermionic degree of freedom in the Standard 
Model (SM), and fermionic partners for Higgs and gauge bosons. So far, the direct search for 
SUSY particles has not been successful. One can only set lower bounds of several hundreds 
of GeV, depending on the particle and the model specifications [3111]. To lift the masses 
of the SUSY partners from the corresponding SM values, soft SUSY-breaking terms are 
introduced [I]. The most general flavor structure of the soft SUSY-breaking sector with 
flavor non-diagonal terms would induce large flavor- changing neutral-currents, contradicting 
the experimental results [5]. Attempts to avoid this kind of problem include flavor-blind 
SUSY-breaking scenarios, like minimal Supergravity or gauge-mediated SUSY-breaking. In 
these scenarios, the sfermion-mass matrices are flavor diagonal in the same basis as the 
quark matrices at the SUSY-breaking scale. However, a certain amount of flavor mixing is 
generated due to the renormalization-group evolution from the SUSY-breaking scale down to 
the electroweak scale. In a more agnostic approach all flavor- violating terms are introduced as 
free parameters, and each model point, i.e. each combination of flavor- violating parameters, 
is tested against experimental data. 

Similarly strong bounds exist for the MSSM Higgs sector from the non-observation of 
Higgs bosons at LEP [6l[7], the Tevatron [8] and most recently from LHC searches [9]. The 
LHC has good prospects to discover at least one neutral Higgs boson over the full MSSM 
parameter space. A precision on the mass of a SM-like Higgs boson of ~ 200 MeV is ex- 
pected [T0j,[Tl] . At the ILC a determination of the Higgs boson properties (within the kine- 
matic reach) will be possible, and an accuracy on the mass could reach the 50 MeV level [T21 
[T3j . The interplay of the LHC and the ILC in the neutral MSSM Higgs sector is discussed 
in Ref. [H]. 

For the MSSvQ the status of higher-order corrections to the masses and mixing angles 
in the neutral Higgs sector is quite advanced. The full one-loop and potentially all leading 
two-loop corrections are known, see Ref. [I7J for a review. Most recently leading three-loop 
corrections became available [T8|[T9"]. 

However, the impact of non-minimal flavor violation (NMFV) on the MSSM Higgs-boson 
masses and mixing angles, entering already at the one-loop level, has not been explored very 
deeply so far. A one-loop calculation taking into account the LL-mixing between the third 
and second generation of scalar up- type quarks has been performed in Ref. [20]. A full 
one-loop calculation of the Higgs-boson self-energies including all NMFV mixing terms had 
been implemented into the Fortran code FeynHiggs [T6l[T7ll2T| l22] . however no cross checks 
or numerical evaluations analyzing the effects of the new mixing terms were performed. 

1 We concentrate here on the case with real parameters. For complex parameters see Refs. [15j[T6] and 
references therein. 
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Possible effects from NMFV on Higgs boson decays were investigated in j23H25j . Within a 
similar context of NMFV there have been also studied some effects of scharm-stop flavor 
mixing in top-quark FCNC processes [2B] and charged Higgs processes [27] as well as the 
implications for LHC [27J|28]. Some previous studies on the induced radiative corrections 
on the Higgs mass from scharm-stop flavor mixing have also been performed in [26l[28], but 
any effects from mixing involving the first generation of scalar quarks have been neglected. 
The numerical estimates in [261 [28] also neglect all the flavor mixings in the scalar down- 
type sector, except for those of LL-type that are induced from the scalar up-type sector by 
SU(2) invariance. In [28] they also consider one example with a particular numerical value 
of non-vanishing sl — mixing. 

We study in this paper the consequences from NMFV for the MSSM Higgs-boson spec- 
trum, where our results are obtained in full generality, i.e. all generations in the scalar up- 
and down-type quark sector are included in our analytical results. In the numerical analysis 
we focus particularly on the flavor mixing between the second and third generations which is 
the relevant one in B physics. Our estimates include all type of flavor mixings, LL, LR, RL^ 
and RR. We devote special attention to the LR/RL sector. These kind of mixing effects are 
expected to be sizable, since they enter the off-diagonal A parameters, which appear directly 
in the coupling of the Higgs bosons to scalar quarks. 

Concerning the constraints from flavor observables we take into account the most up-to- 
date evaluations in the NMFV MSSM for BR(B -»■ X s j% BR(B S -»■ and AM Bs , 
based on the BPHYSICS subroutine included in the SuFla code |31j . For the evaluation of 
AMb s we have incorporated into this subroutine additional contributions from the one-loop 
gluino boxes [52] which are known to be very relevant in the context of NMFV scenarios [551 - 
[55]. 

In the first step of the analysis we scan over the NMFV parameters, and contrast them 
with the experimental bounds on BK(B — >■ X a j), BK(B S — > and AMg s . In the 

second step we analyze the one-loop contributions of NMFV to the MSSM Higgs boson 
masses, focusing on the parameter space still allowed by the experimental flavor constraints. 
In this way the full possible impact of NMFV in the MSSM on the Higgs sector can be 
explored. The paper is organized as follows: in Sect. [2] we introduce our notation for the 
NMFV MSSM and define certain benchmark scenarios that are used for the subsequent 
analysis. Our implementation and new results on i?-physics observables is given in Sect. [5J, 
where we also analyze in detail which combination of NMFV parameters are still allowed by 
current experimental constraints. The calculation of the corrections to Higgs boson masses 
in the NMFV MSSM is presented in Sect. [U The numerical analysis of the impact of the 
one-loop Higgs mass corrections is given in Sect. [5j Our conclusions can be found in Sect. [61 

2 SUSY scenarios with Non- Minimal Flavor Violation 

We work in SUSY scenarios with the same particle content as the MSSM, but with general 
flavor mixing hypothesis in the squark sector. Within these SUSY-NMFV scenarios, besides 
the usual flavor violation originated by the CKM matrix of the quark sector, the general flavor 

2 Subleading NLO MSSM corrections were evaluated in [29 1 130 ] . However, their effect on our evaluations 
would be minor. 



2 



mixing in the squark mass matrices additionally generates flavor violation from the squark 
sector. These squark flavor mixings are usually described in terms of a set of dimensionless 
parameters 5* Y (X, Y = L,R; i, j = t,c,u or b, s, d), which for simplicity in the computations 
are frequently considered within the Mass Insertion Approximation (MIA) [36] . We will not 
use here this approximation, but instead we will solve exactly the diagonalization of the 
squark mass matrices. 

In this section we summarize the main features of the squark flavor mixing within the 
SUSY-NMFV scenarios and set the notation. 
The relevant MSSM superpotential terms are: 



W 



H%QPY U U - Q l3 Y d D + (iH?H% 



where the involved superfields are: Q, containing the quark (ul di) T and squark (ul di) T 
SU(2) doublets; U, containing the quark {ur) c and squark u* R SU{2) singlets; D, containing 
the quark (dn) c and squark d* R 577(2) singlets; and Hi j2 containing the Higgs bosons SU(2) 
doublets, Hi = (H? ) T and T-L 2 = {H 2 H 2 ) T i an d their SUSY partners. f c denotes here 
the particle-antiparticle conjugate of fermion /, and /* denotes the complex conjugate of 
sfermion /. We follow the notation of [37], but with the the convention e 12 = —1. The 
generation indices in the superfields, Qi,Ui, £>i, quarks qi, squarks q iy (i = 1,2,3), and 
Yukawa coupling 3x3 matrices, Y£, Y&, = 1, 2, 3), have been omitted above for brevity. 

Usually one starts rotating the quark fields from the SU (2) (interaction) eigenstate basis, 
<?Li?,5 to the mass (physical), ^ h ^ s eigenstate basis by unitary transformations, V^ R : 



U L,R 
„phys 



\ 



phys 



V L,R 



U L,R 
„int 
C L,R 
j.int 
L,R 



( 



\ 



dt y R 

„phys 
b L,R 
z,phys 
U L,R 



\ 



7 



yd 



Jint 

a L,R 

c int 

b L,R 

tint 

°L,R 



(2) 



such that the fermion mass matrices in the physical basis are diagonal: 



V£Y U *V R 



V T d Y 



R 



diag(2/ u , y c , Vt) 
diag(y d , y„, y b ) 



( m u m c m t 
diag — , — , — 

V v 2 v 2 v 2 



diag 



(3) 
(4) 



m d m s m b 

Vi ' V i ' Ui 

where V\ = (H\) and v 2 = (H 2 ) are the vacuum expectation values of the neutral Higgs 
fields. The CKM matrix, which is responsible for the flavor violation in the quark sector, is 
given as usual as, 

Vckm = VZV?. (5) 
A simultaneous (parallel) rotation of the squarks with the same above unitary matrices as 
their corresponding quark partners leads to the so-called Super-CKM basis. In the NMFV 
scenarios, contrary to the MFV ones, the Super-CKM basis does not coincide with the 
physical squark basis, i.e, their corresponding squark mass matrices are not yet diagonal. 
More concretely, rotating the squarks from the interaction basis, q 1 ^ to the Super-CKM 
basis, q L R , by 




~ int 



yu 



U L,R 
~int 

l. L ' R 
flat 
L,R 




L,R 



Jint 

a L,R 

~int 

t L > R 

Lint 

°L,R 



(6) 
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one gets the soft-SUSY-breaking Lagrangian transformed from: 

Q^A\Uf^U 2 - QfA%Vf*U 1 + h.c. 

to: 



(7) 



£ S oft = -U Ri m% R .fl Rj - Vyn% Rij V Rj - Ul#n\ hi U Li - V* Li m% Lij V Lj (8) 
- KuAijUxjl-L® - V Li {VcKM) kiAljU^T-i^ ~ ^Li(Vc KM )ikA kj T> Rj Hi + V , Li Af j V* Rj 'Hi + h.c. 

where we have used calligraphic capital letters for the squark fields with generation indexes, 

;"/int _ ~int ~int Tint. -y^int _ Tint ~int lint Aint _ /~int Jint\T /~int ~int\T /Tint %mt\T . /q\ 
U l,2,3 — U R i c R i l R i U \,2,3 ~ a R i s R i°R i ^1,2,3 — \ U L a L ) i\ c L S L ) i \ l L °L ) ) l y J 

Wli,2,3 = ml, c l , t L ; £^1,2,3 = d L , §l, b L ; U R1:2j3 = Ur, c r , t R ; T> Rh2:3 = d R , s R: b R ; (10) 
and (q = u, d) 

A" = VlA«Vg,ml R = VlmlVl\m\ R = V^V^m 2 ^ = y?m|y* = V^mJV?. 

(11) 

The usual procedure to introduce general flavor mixing in the squark sector is to include the 
non-diagonality in flavor space at this stage, namely, in the Super-CKM basis. Thus, one 
usually writes the 6x6 non-diagonal mass matrices, M.\ and M.% referred to the Super- 
CKM basis, being ordered respectively as (ul, cl, tL, u R , c R , t R ) and (o^, SL,bi,d R , s R ,b R ), 
and write them in terms of left- and right-handed blocks M~ AB (q = u,d, A, B = L,R), 
which are non-diagonal 3x3 matrices, 

M}=(S L « = «U (12) 

where: 



V Mf LR M~ RR 



M*LLij = m l Lij + « + ( T s ~ Q«sin 2 MMf cos2/3) <5y, 
Mi RRlJ =ml R .. + (m 2 Ui + Q u sm 2 6 w M 2 z cos2(3) 
Mi LRij = (n° 2 ) A% - m Uif i cot /3 Sij, , 
M iLLij = m kij + « + ( T " - Qd ^ 0w)M 2 z cos 2/3) S ij} 
M SRRij = m l Rij + ( m l + sin 2 6 W M Z cos 2/3) 

M l LRij =W) A%-m Aii i tan/3%. (13) 

with, i,j = 1,2,3, Q u = 2/3, Q d = -1/3, T 3 U = 1/2 and T 3 d = -1/2. 6 W is the weak angle, 
M z is the Z gauge boson mass, and (m Ul ,m U2 ,m U3 ) = (m u ,m c ,m t ), (m dl ,m d2 ,m d3 ) = 
(m d , m s , rrtfe). It should be noted that the non-diagonality in flavor comes from the values of 
m 2 - T . ., m 2 ~ T . ., rrA . ., m\ . ., Af 4 and Af, for % ^ j. 
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The next step is to rotate the squark states from the Super-CKM basis, c[l,r, to the phys- 
ical basis, q phys . If we set the order in the Super-CKM basis as above, (ul, cl, ti, ur, cr, tn) 
and (d,L, §l, Sl, (Ir, §r, 1>r), and in the physical basis as u li 6 and di r _ 6 , respectively, these last 
rotations are given by two 6x6 matrices, R u and R d , 



( ui \ 

U i 

V % / 



R u 



( u L \ 

h 

UR 

cr 
V ** ) 



( di \ 

d-2 
d 3 
d>4 
d 5 
\ d 6 ) 



R a 



( d L \ 

sl 
h 
d~R 
Sr 
\~b R J 



(14) 



yielding the diagonal mass-squared matrices as follows, 



diagjm^ , m\ 2 , m? 3 , m\ , m? 5 , m u& 



diagjmj , 



m 2 -, , m j , m j , m j , m , 

d3 ' d 5 ' da 



"d 2 > 



2 } 



St 



R u Ml R 



R d M\ R d ^ . 

a 



(15) 
(16) 



The corresponding Feynman rules in the physical basis for the vertices needed for our 
applications, i.e. the interaction of one and two Higgs or gauge bosons with two squarks, 
can be found in the Appendix A. This new set of generalized vertices had been implemented 
into the program packages FeynArts/ FormCalc [381139] extending the previous MSSM model 
file [10] • The extended FeynArts version was used for the evaluation of the Feynman diagrams 
along this paper to obtain the general analytical results. 

In the numerical part of the present study we will restrict ourselves to the case where 
there are flavor mixings exclusively between the second and third squark generation. These 
mixings are known to produce the largest flavor violation effects in B meson physics since 
their size are usually governed by the third generation quark masses. On the other hand, 
and in order to reduce further the number of independent parameters, we will focus in the 
following analysis on constrained SUSY scenarios, where the soft mass parameters fulfill 
universality hypothesis at the gauge unification (GUT) scale. Concretely, we will work 
with the so-called Constrained MSSM (CMSSM) and Non Universal Higgs Mass (NUHM) 
scenarios, which are defined by (see [12] and references therein), 



CMSSM : 
NUHM : 



m , mi/2, A , sign(/i), tan/3 
m , mi/2, A , sign(/x), tan 0, m Hl ,m H2 



(17) 



where, A is the universal trilinear coupling, m , mi/2, m^, rnn 2 , are the universal scalar 
mass, gaugino mass, and Hi and H 2 Higgs masses, respectively, at the GUT scale, sign(/x) is 
the sign of the /x parameter and again tan/3 = v 2 jv\. The soft Higgs masses in the NUHM 
are usually parametrized as rn 2 Hi2 = (1 + Si^m^, such that by taking 5i 2 = one moves 
from the NUHM to the CMSSM.' 



It should be noted that the condition of universal squark soft masses, m 2 - 



i i 
m- . . = m% . . 

Dm Drij 



2 



m^Sij, is fulfilled just at the GUT scale. When running these soft mass 
matrices from the GUT scale down to the relevant low energy, they will generically turn 
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non-diagonal in flavor. However, in MFV scenarios the non-diagonal terms are exclusively 
generated in this running by off-diagonal terms in the Vckmj and therefore they can be 
safely neglected at low energy. Contrary, in NMFV scenarios, the universal hypothesis in 
these squark mass matrices is by definition not fulfilled at low energies. 

Our final settings for the numerical evaluation of the squark flavor mixings in NMFV 
scenarios are fixed (after RGE running) at low energy as follows, 



/ m 



m 



2 

Ulu 










m fJ 

r-f J UL22 

^23 m L22 m L33 



o 



S 23 Lr ^U L22 m U L33 

m L 3 



v 2 A u 









t m U R 22 m U L33 



U R11 







m t A t 




U R 22 

c * U R 22 Ur 33 



5% R mf T m T ~ T 

Ct U R 22 U R33 



m 



u R33 



\ 
J 



Vl A a 



in 



Of 



( 



\ 



rn 






5 fb L ^b R 22 m D L33 







2 










2 


D R 22 





S^m f, m 

SO U R 22 



S so R ^D L 22 m D R33 

m h A h 




°sb 



L> R 22 L > R33 

m 



D R33 



(19) 



(20) 



[21] 



(22) 



(23) 



It is worth mentioning that the relation between the two soft squark mass matrices in the 
'Left' sector ( 12T|) is due to SU(2) gauge invariance. Eq. [21] can be derived from the two 
last relations of eq. [TTJ This dependence between the non-diagonal terms of these squark 
mass matrices is the reason why is introduced 5%$ instead of two independent deltas S^ t L 
and 5g b L . To get the needed running of the soft parameters from the GUT scale down to 
low energy, that we set here 1 TeV, we solve numerically the one-loop RGEs with the code 
SPHENO jl3]. The diagonalization of all the mass matrices is performed with the program 
FeynHiggs HSHI3EI1I22] • 

In CMSSM and other SUSY-GUT scenarios the flavor changing deltas go (in the leading 

logarithmic approximation) as 5^ — ^"^V 4 ^ (Y q ^Y q )23 log( ^^ ) (m 2 is usually taken 
as the geometric mean of the involved flavor diagonal squared squark mass matrix entries, 
see eq. El]), whereas the LR, RL and RR deltas are suppressed instead by small mass ratios, 



{iriqAp) 



and 



K) 



respectively. Furthermore, in these scenarios the mixing involving the 



first generation squarks is naturally suppressed by the smallness of the corresponding Yukawa 
couplings. In order to keep the number of free parameters manageable, this motivated our 
above choice of neglecting in the numerical analysis the mixing of the first generation squarks. 
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However, we will not assume any explicit hierarchy in the various mixing terms between the 
second and third generation. 

It should be noted that in the 'Left-Right' sector we have normalized the trilinear cou- 
plings at low energies as A\a = VqA^ (with A 33 = A t and A 33 = A b ) and we have neglected 
the Ai couplings of the first and second generations. Finally, it should be noted that the 
dimensionless parameters 5^ Y defining the non-diagonal entries in flavor space (i ^ j) are 
normalized respect the geometric mean of the corresponding diagonal squared soft masses. 
For instance, 

Sj* = (v 2 A u )23/(mu L22 ™u R Ji *ct L = (v2A u ) 32 /(m &R2 m &L33 ), etc. (24) 

For definiteness and simplicity, in the present work we will perform our estimates in 
specific constrained SUSY scenarios, CMSSM and NUHM, whose input parameters m , 7711/2, 
A , tan/3, sign(/i), 5 12 , are summarized in tabled, and supplemented with bf Y as described 
above. Regarding CMSSM, we have chosen six SPS benchmark points [2], SPSla, SPSlb, 
SPS2, SPS3, SPS4, and SPS5 and one more point with very heavy spectra (VHeavyS). It 
should be noted that several of these SPSX points are already in conflict with recent LHC 
data [Hill], but we have chosen them here as reference points to study the effects of SUSY 
on the Higgs mass corrections, since they have been studied at length in the literature. 
At present, a heavier SUSY spectrum, as for instance our point VHeavyS is certainly more 
realistic and compatible with LHC data. In general an analysis of LHC data including NMFV 
effects in the squark sector would be very desirable. Regarding NUHM, we have chosen a 
point with heavy SUSY spectra and light Higgs sector ( Heavy SLightH) and a point (BFP) 
that has been proven in [15] to give the best fit to the set of low energy observables. For later 
reference, needed in our posterior analysis of the Higgs mass corrections, we also include in 
the table the corresponding MSSM Higgs masses, computed with FeynHiggs [T6 | [T7 ][2Tl[22] 
and with all flavor changing deltas switched off, i.e., 5^ Y = 0. 



3 Constraints on Non- Minimal Flavor Violating SUSY 
scenarios from _B-Physics 

In this section we analyze the constraints on Non-Minimal Flavor Violating SUSY scenarios 
from Z?-Physics. Since we are mainly interested in the phenomenological consequences of the 
flavor mixing between the third and second generations we will focusj on the following three B 
meson observables: 1) Branching ratio of the B radiative decay BR(I? — > X s j), 2) Branching 
ratio of the B s muonic decay BK(B S — > and 3) B s — B s mass difference AMb s - 

3 We adopt here the definition in terms of the GUT-scale input parameters, while the original definition 
in [33] was based on the weak-scale parameters. 

4 We have checked that electroweak precision observables, where NMFV effects enter, for instance, via 
[20], do not lead to relevant additional constraints on the allowed parameter space. Our results on this 
constraint are in agreement with Ref. [26]. 
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points 


mi/j 


m 


A 


tan {3 


6i 


§2 


m h 


m H 


M A 


m H ± 


SPSla 


250 


100 


-100 


10 








112 


394 


394 


402 


SPSlb 


400 


200 





30 








116 


526 


526 


532 


SPS2 


300 


1450 





10 








115 


1443 


1443 


1445 


SPS3 


400 


90 





10 








115 


573 


572 


578 


SPS4 


300 


400 





50 








114 


404 


404 


414 


SPS5 


300 


150 


-1000 


5 








111 


694 


694 


698 


VHeavyS 


800 


800 


-800 


5 








120 


1524 


1524 


1526 


HeavySLightH 


600 


600 





5 


-1.86 


+ 1.86 


114 


223 


219 


233 


BFP 


530 


110 


-370 


27 


-84.7 


-84.7 


120 


507 


507 


514 



Table 1: Values of m 1( / 2 , m , A , tan/3, 8i, 5 2 and Higgs boson masses m h , m H , M4 and 
m,H± for the points considered in the analysis. All parameters with mass dimension are in 
GeV, and sign(/i) > for all points. 



Another B observable of interest in the present context is BR(5 — > X s l + l~). However, we 
have not included this in our study, because the predicted rates in NMFV-SUSY scenarios 
for this observable are closely correlated with those from BK(B — > X s , ~f) due to the dipole 
operators dominance in the photon-penguin diagrams mediating BR(I? — > X s l + l~) decays. It 
implies that the restrictions on the flavor mixing 5* Y parameters from BR(i? — > X s l + l~)are 
also expected to be correlated with those from the radiative decays. 

The summary of the relevant features for our analysis of these three B meson observables 
is given in the following. 



3.1 BR(£? X s 7) 

The relevant effective Hamiltonian for this decay is given in terms of the Wilson coefficients 
Ci and operators 0, by: 

AC 8 

Kfi = -^f^KM^CKM + C i°^ ( 25 ) 

Where the primed operators can be obtained from the unprimed ones by replacing L -h- R. 
The complete list of operators can be found, for instance, in j4T] . In the context of SUSY 
scenarios with the MSSM particle content and MFV, only two of these operators get relevant 
contributions, the so-called photonic dipole operator O7 and gluonic dipole operator 0% given, 
respectively, by: 

Oi = -^- 2 m b {s L a^b R )F^ (26) 



Os = ^m b (s L a^T% R ) G^. (27) 

We have omitted the color indices here for brevity. Within NMFV also 0' 7 8 have to be taken 
into account: 



0' 7 = ^m b (s R a^b L )F^, (2* 



O's = ^m b (s R a^T% L ) G%. (29) 
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The Wilson coefficients at the SUSY scale are obtained as usual by the matching proce- 
dure of the proper matrix element being computed from the previous effective Hamiltonian 
to the corresponding matrix elements being computed from the SUSY model operating at 
that SUSY scale, the NMFV-MSSM in our case. These Wilson coefficients encode, there- 
fore, the contributions to BR(B — » X s ^) from the loops of the SUSY and Higgs sectors of 
the MSSM. The effects from squark flavor mixings that are parametrized by the 5* Y , are 
included in this observable via the squark physical masses and rotation matrices, given in 
the previous section, that appear in the computation of the matrix element at the one loop 
level and, therefore, are also encoded in the Wilson coefficients. The explicit expressions for 
these coefficients in the MSSM, in terms of the physical basis, can be found, for instance, 
in refs. |^6T - H8] . We have included in our analysis the most relevant loop contributions to 
the Wilson coefficients, concretely: 1) loops with Higgs bosons, 2) loops with charginos and 
3) loops with gluinos. It should be noted that, at one loop order, the gluino loops do not 
contribute in MFV scenarios, but they are very relevant (dominant in many cases) in the 
present NMFV scenarios. 

Once the Wilson coefficients are known at the SUSY scale, they are evolved with the 
proper Renormalization Group Equations (RGEs) down to the proper low-energy scale. As 
a consequence of this running the previous operators mix and the corresponding Wilson 
coefficients, Cj$ get involved in the computation of the B — > X s ^ decay rate. The RGE- 
running is done in two steps: The first one is from the SUSY scale down to the electroweak 
scale, and the second one is from this electroweak scale down to the B-physics scale. For 
the first step, we use the LO-RGEs for the relevant Wilson coefficients as in [48] and fix six 
active quark flavors in this running. For the second running we use the NLO-RGEs as in [49] 
and fix, correspondingly, five active quark flavors. For the charged Higgs sector we use the 
NLO formulas for the Wilson coefficients as in [50] . 

The resummation of scalar induced large tan/3 effects is performed, as usual, by the 
effective Lagrangian approach that parametrizes the one-loop generated effective couplings 
between the "H 2 Higgs doublet and the down type quarks in softly broken SUSY models [5T] . 
A necessary condition to take into account all tan /3-enhanced terms in flavor changing am- 
plitudes is the diagonalization of the down-type quark mass matrix in the presence of these 
effective couplings [52TI54"] . A summary of this effective Lagrangian formalism for the resum- 
mation of large tan j3 effects in the three B observables of our interest, within the context of 
MFV scenarios, can be found in [55]. We follow here the treatment of [56] where the resum- 
mation of large tan /3 effects via effective Lagrangians is generalized to the case where the 
effective d^dj^H^ coupling contains also non-minimal sources of flavor mixing. It should be 
noted that the most relevant scalar induced large tan f3 effects for the present work are those 
generated by one-loop diagrams with gluino-sbottom and chargino-stop inside the loops. 
The large tan /3 resummation effects and the relevance of other chirally enhanced corrections 
to FCNC processes within the NMFV context have recently been studied exhaustively also 
in [57J|5B] (previous studies can be found, for instance, in Refs. [59~HB"T] ). 

The total branching ratio for this decay is finally estimated by adding the new contri- 
butions from the SUSY and Higgs sectors to the SM rate. More specifically, we use eq.(42) 
of [49] for the estimate of BK(B — > A s 7) in terms of the ratios of the Wilson coefficients 
CV 5 8 and C' 7 8 (including all the mentioned new contributions) divided by the corresponding 
Cff in the SM. 
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For the numerical estimates of BR(£> — > X s ^) we use the FORTRAN subroutine BPHYSICS 
(modified as to include the contributions from C 7 8 which were not included in its original ver- 
sion) included in the SuFla code, that incorporates all the above mentioned ingredients [3"T] . 

Finally, for completeness, we include below the experimental measurement of this ob- 
servable [5l l62pl . and its prediction within the SM [63] : 

BR(B X s7 ) cxp = (3.55 ± 0.26) x 10" 4 (30) 
BR(B X s7 ) SM = (3.15 ± 0.23) x 10" 4 (31) 



3.2 BR(B S -> 



The relevant effective Hamiltonian for this process is [63] [65] : 



(32) 



where the operators are given by: 

Oio = (s^PLb) {Plulsv) 
O s = m b (sPpb) (pp) , 
P = m b (sP R b) (/275A«) , 



0' s = m s (sP L b) (p/S) , 
O'p = m s (sP L b) (prfon) . 



(33) 



We have again omitted the color indices here for brevity. 

In this case, the RG running is straightforward since the anomalous dimensions of the 
above involved operators are zero, and the prediction for the decay rate is simply expressed 
by: 



s~~i'2 2 2 £2 

Tyr> I D , ,,+ ,,-\ _ Lr F a m B s JB s T B s ^ rtsif T/t6 | 2 f, TT~ 2 

BR(# s ^/i p )- I^ckm^ckmI y 1 " 4 ^ 



x 



[(l-4ml)\F s \ 2 + \F P + 2mlF l0 \ 2 } 
where m M = m VL jmB s and the Fi are given by 

~C s ,pm b - C' SP m s 



(34) 



S,P 



m b + m s 



'10 



C 
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Within the SM the most relevant operator is 0\q as the Higgs mediated contributions to 
Os,p can be safely neglected. It should be noted that the contribution from Oio to the decay 
rate is helicity suppressed by a factor of since the B s meson has spin zero. In contrast, 
in SUSY scenarios the scalar and pseudo-scalar operators, Os,p, can be very important, 
particularly at large tan/3 > 30 where the contributions to Cs and Cp from neutral Higgs 
penguin diagrams can become large and dominate the branching ratio, because in this case 
the branching ratio grows with tan/3 as tan 6 /3. The studies in the literature of these MSSM 
Higgs-penguin contributions to BR(B S — >■ /i + /i~) have focused on both MFV [5lJ[66l[67] an d 
NMFV scenarios (33JE51EI1EI] • In both cases the rates for BR(_B S — > fi + fi~) at large tan/3 



3 We have added the various contributions to the experimental error in quadrature. 
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can be enhanced by a few orders of magnitude compared with the prediction in the SM, 
therefore providing an optimal window for SUSY signals. 

In the present context of SUSY-NMFV, with no preference for large tan /3 values, there 
are in general three types of one-loop diagrams that contribute to the previous Cj Wilson 
coefficients for this B s — > fi + fi~ decay: 1) Box diagrams, 2) Z-penguin diagrams and 3) 
neutral Higgs boson if-penguin diagrams, where H denotes the three neutral MSSM Higgs 
bosons. In our numerical estimates we have included what are known to be the dominant 
contributions to these three types of diagrams [M]: chargino contributions to box and Z- 
penguin diagrams and chargino and gluino contributions to if-penguin diagrams. 

Regarding the resummation of large tan f3 effects we have followed the same effective 
Lagrangian formalism as explained in the previous case of B — > X s j. In the case of contri- 
butions from if -penguin diagrams to B s — > fi + fi~ this resummation is very relevant and it 
has been generalized to NMFV-SUSY scenarios in |56j . 

For the numerical estimates we use again the BPHYSICS subroutine included in the 
SuFla code [31] which incorporates all the ingredients that we have pointed out above. 

Finally, for completeness, we include below the present experimental upper bound for 
this observable [68], and the prediction within the SM |69j : 

BR(B S ^ + /i") CX p < 1.1 x 1(T 8 (95% CL) (35) 
BR(B S /i + /OsM = (3.6 ± 0.4) x 10" 9 (36) 



3.3 AM Bs 

The relevant effective Hamiltonian for B s — B s mixing and, hence, for the B s /B s mass 
difference AMb s is: 

n cS = ^M 2 W (V$£ U V& U ) 2 CiOi. (37) 

i 

In the SM the most relevant operator is: 

vll = (6 ViS «)^ 7 M P ^). ( 38) 

Where we have now written explicitly the color indices. In scenarios beyond the SM, as the 
present NMFV MSSM, other operators are also relevant: 

lr = (p« 1 ^p L s a )(W^P R s p ), 0\ R = (b a P L s a )(¥P R s' 3 ), (39) 

sll = (6«p LS «)(^p^) ) O s 2 LL = {b a a^P L s a )(¥a^P L s p ), (40) 

and the corresponding operators VRR and Of RR that can be obtained by replacing P^ -h- Pr 
in fl38|) and (T40]) . The mass difference AMb s is then evaluated by taking the matrix element 

AM Bs =2\(B s \n cS \B s }\, (41) 
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where (B s \TL e ff\B s ) is given by 

(B s \H eS \B s ) =^M 2 w m B J 2 Ba (C^ckm) 2 E W ■ (42) 

j 

Here mg s is the _B S meson mass, and fs s is the B s decay constant. The coefficients Pi contain 
the effects due to RG running between the electroweak scale fiw and m& as well as the relevant 
hadronic matrix element. We use the coefficients Pi from the lattice calculation [70] : 

pVLL =0 .73 7 pLR = _ i.97 ) pLR =2.50, Pf LL = - 1.02, P 2 5LL = - 1.97. (43) 

The coefficients Pj , etc., may be obtained from those above by simply exchanging L ^ R. 

In the present context of SUSY-NMFV, again with no preference for large tan /3 values, 
and besides the SM loop contributions, there are in general three types of one-loop diagrams 
that contribute to the previous Cj Wilson coefficients for B s — B s mixing: 1) Box diagrams, 
2) Z-penguin diagrams and 3) double Higgs-penguin diagrams. In our numerical estimates 
we have included what are known to be the dominant contributions to these three types of 
diagrams in scenarios with non-minimal flavor violation (for a review see, for instance, [33]): 
gluino contributions to box diagrams, chargino contributions to box and Z-penguin diagrams, 
and chargino and gluino contributions to double i7-penguin diagrams. As in the previous 
observables, the total prediction for AMb 3 includes, of course, the SM contribution. 

Regarding the resummation of large tan j3 effects we have followed again the effective 
Lagrangian formalism, generalized to NMFV-SUSY scenarios [SB], as in the previous cases 
of B — > Xg'y and B s — > /i 4 pT . It should be noted that, in the case of AMg s , the resummation 
of large tan /3 effects is very relevant for the double if-penguin contributions, which grow 
very fast with tan/3. 

For the numerical estimates we have modified the BPHYSICS subroutine included in the 
SuFla code [3T] which incorporates all the ingredients that we have pointed out above, except 
the contributions from gluino boxes. These contributions are known to be very important for 
B s — B s mixing in SUSY scenarios with non-minimal flavor violation [33H35] and therefore 
they must be included into our analysis of AMb s - Concretely, we have incorporated them 
into the BPHYSICS subroutine by adding the full one-loop formulas for the gluino boxes 
of [32] to the other above quoted contributions that were already included in the code. In 
order to illustrate the relevance of these gluino contributions in our analysis of AMb s , we 
show in Fig{T] each separate contribution as a function of tan /3 in a particular example with 
$23 = $rf) R = 0.1, that we have chosen for comparison with [33]. The other flavor changing 
deltas are set to zero, and the other relevant MSSM parameters are set to rriq = 500 GeV, 
A t = —rrig, m- g = y/2rriq, \x = nig, and tua = 300 GeV, as in Fig. 24 of [33]. We clearly 
see in FigJT] that it is just in the very large tan f3 region where double Higgs- penguins 
dominate. For moderate and low tan (3 values, tan f3 < 20 (which is a relevant region for 
our numerical analysis, see below) the gluino boxes fully dominates the SUSY corrections 
to AMb s and compete with the SM contributions. Our numerical estimate in this plot is in 
complete agreement with the results in [33] (see, in particular, Fig. 24 of this reference) which 
analyzed and compared in full detail these corrections. Finally, for completeness, we include 
below the experimental measurement of this observable [5_J§, and its prediction within the 

6 We have again added the various contributions to the experimental error in quadrature. 
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Figure 1: Relevant contributions to AM# s in NMFV-SUSY scenarios as a function of tan/3. 
They include: SM, Double Higgs penguins, gluino boxes and chargino boxes. The total 
prediction for AMg s should be understood here as AMb s = |Total|. The parameters are set 
to 623 = 5 RR = 0.1,771, = 500 GeV, A t = -m,-, rrig = y/2rrig, [i = m,-, and niA = 300 GeV. 
The other flavor changing deltas are set to zero. 



SM (using NLO expression of [7T] and error estimate of [72]): 

AM B,cx P = (H7.0 ± 0.8) x 10^ 10 MeV , (44) 
AM BaSM = (117.1^) x 10- 10 MeV . (45) 

3.4 Numerical results on B observables 

In the following of this section we present our numerical results for the three B observables 
in the NMFV-SUSY scenarios and a discussion on the allowed values for the flavor changing 
deltas, &f Y - 

The predictions for BR(£> — > X s j), BR(B S — >■ )and AMg s versus the various bf Y , 
for the six selected SPS points, are displayed respectively in Figs. El [3] and HI For this 
analysis, we have assumed that just one at a time of these deltas is not vanishing. Results 
for two non-vanishing deltas will be shown later. The following 7 flavor changing deltas are 
considered: <5^ 3 L , 5^ R , 5^ b R , 5 RL , 5 R b L , 5 RR and 5 RR ; and we have explored delta values within 
the interval — 1 < bf Y < 1. In all plots, the predictions for 5* Y = represent our estimate 
of the corresponding observable in the MFV case. This will allow us to compare easily the 
results in the two kind of scenarios, NMFV and MFV. It should be noted also, that some of 
the predicted lines in these plots do not expand along the full interval —1 < bf Y < 1, and 
they are restricted to a smaller interval; for some regions of the parameter space a too large 
delta value can generate very large corrections to any of the masses, and the mass squared 
turns negative. These problematic points are consequently not shown in our plots. 
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We have also included in the right vertical axis of these figures, for comparison, the 
respective SM prediction in ( 1HT1) . (1561) . and (J4"5]h The error bars displayed are the corre- 
sponding SM uncertainties as explained below. The shadowed horizontal bands in the case 
of BR(I? — > Xg 1 -)) and AMg s are their corresponding experimental measurements in (1301 ). 
and (JHJ), expanded with 3cr exp errors. In the case of BR(£> S — > the shadowed area 

corresponds to the allowed region by the upper bound in (1351) . 

The main conclusions extracted from these figures for the three B observables are sum- 
marized as follows: 

• BR(5 X sl ): 

- Sensitivity to the various deltas: 

We find strong sensitivity to 5g b R , 5 R b L , 5^, S^ R and S^ R , in all the studied points, 
for both high and low tan /3 values. The order found from the highest to the 
lowest sensitivity is, generically: 1) Sf R and 5 R b L the largest, 2) 5^ the next, 3) 
S^ R and S RR the next to next, and 4) slight sensitivity to 5 RR and 5 RL . 

- Comparing the predictions with the experimental data: 

If we focus on the five most relevant deltas, 5^ R , 5 R b L , 5%% , S RR and 5^ R , we see 
clearly that tiny deviations from zero (i.e., deviations from MFV) in these deltas, 
and specially in the first three, produce sizeable shifts in BK(B — > X s j), and in 
many cases take it out the experimental allowed band. Therefore, it is obvious 
from these plots that BR(_B — > X s j) sets stringent bounds on the deltas (when 
varying just one delta), which are particularly tight on 5^ R , 5 R b L , 5%%, 5 R b R , and 
S^ R , indeed for all the studied SPS points. There are just two exceptions, where 
the predicted central values of BR(I? — > X s '-f) are already outside the experimental 
band in the MFV case (all deltas set to zero), and assuming one of these three most 
relevant deltas to be non- vanishing, the prediction moves inside the experimental 
band. This happens, for instance, in the points SPS4 and SPSlb that have the 
largest tan (3 values of 50 and 30 respectively. Interestingly, it means that some 
points of the CMSSM, particularly those with large tan (3 values, that would have 
been excluded in the context of MFV, can be recovered as compatible with data 
within NMFV-SUSY scenarios. 

- Intervals of 5* Y allowed by data: 

In order to conclude on the allowed delta intervals we have assumed that our pre- 
dictions of BR(i? — > Xs'-f) within SUSY scenarios have a somewhat larger theo- 
retical error A thco (BR(£? X s7 )) than the SM prediction A t s ^°(BR( J B Ayy)) 
given in ( 13T|) . As a very conservative value we use A theo (BR(i? — > X s ~f)) = 
0.69 x 10~ 4 . A given 5^ Y value is then considered to be allowed by data if the cor- 
responding interval, defined by BR(5 — > X^) ± A thco (BR(i? — > X s ^y)), intersects 
with the experimental band. It corresponds to adding linearly the experimental 
uncertainty and the MSSM theoretical uncertainty. In total a predicted ratio in 
the interval 

2.08 x 10^ 4 < BR(B X s7 ) < 5.02 x 10" 4 , (46) 
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is regarded as allowed. Our results for these allowed intervals are summarized 
in table |2J In this table we see again that the less constrained parameters by 
BR(B — > X s ^) are 5 RL , and 5 RR . Therefore, except for the excluded SPS4 case, 
these two deltas can be sizeable, \S^ Y \ larger than 0(0.1), and compatible with 
BR(B X s7 ) data. 

• BR(B S 

- Sensitivity to the various deltas: 

We find significant sensitivity to the NMFV deltas in scenarios with very large 
tan/3, as it is the case of SPS4 and SPSlb. This sensitivity is clearly due to the 
Higgs-mediated contribution that, grows as tan 6 /3. The largest sensitivity is to 
623- In the case of SPS4, there is also significant sensitivity to 5^ R , 5 RR and S^ R . 
In the SPSlb scenario there is also found some sensitivity to Sf b R , 5 RR , 5 RR and 

xLR 
°ct ■ 

- Comparing the predictions with the experimental data: 

Fig. |3] clearly shows that most of the 1 bf Y \ < 1 explored values are allowed by 
experimental data on BR(B S — > It is in the points with very large tan /3, i.e 

SPS4 and SPSlb, where there are some relevant differences between the MFV and 
the NMFV predictions. First, all predictions for MFV scenarios except for SPS4, 
are inside the experimental allowed area. In the case of SPSlb, the comparison 
of the NMFV predictions with data constraints specially 5^, but also S^ R , 5 RR , 
5 R b R and 5^ R . In the case of SPS4 some input non- vanishing values of 5^, 5g R or 
5 RR can place the prediction inside the experimental allowed area. In the case of 
the SPSla and SPS3 scenarios some constraints for 5^ can be found. 

- Intervals of &f Y allowed by data: 

As in the previous observable, we assume here that our predictions for BR(B S — > 
have a slightly larger error as the SM prediction in ( )36|) . where, how- 
ever, the theory uncertainty is very small in comparison with the experimen- 
tal bound. We choose A thco (BR( J B s = 0.12 x 10~ 8 . Then, a given 
bf Y value is allowed by data if the predicted interval, defined by BR(£> S — >■ 
yU + /i _ ) + A thco (BR(£> s — > intersects the experimental area. The upper 
line of the experimental area in this case is set by the 95% CL upper bound given 
in f )35|) . It implies that for a predicted ratio to be allowed it must fulfill: 

BR(B S < 1.22 x 10" 8 . (47) 

The results for the allowed 5* Y intervals are collected in table [2j We conclude from 
this table that, except for scenarios with large tan/3 > 30, like SPS4 and SPSlb, 
the size of the deltas can be sizeable, | larger than 0(0.1), and compatible 
with BK(B S — > data. 

• AM Bs : 
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- Sensitivity to the various deltas: 

Generically, we find strong sensitivity to most of the NMFV deltas in all the 
studied points, including those with large and low tan/3 values. The pattern of 
the AMb s predictions as a function of the various 5^ Y differs substantially for 
each SPS point. This is mainly because in this observable there are two large 
competing contributions, the double Higgs penguins and the gluino boxes, with 
very different behavior with tan /3, as we have seen in Fig. [TJ In the case of SPS4 
with extremely large tan/3 = 50, the high sensitivity to all deltas is evident in 
this figure. In the case of SPS5 with low tan /3 = 5, there is important sensitivity 
to all deltas, except S RR , S^ R and 5 RL . Generically, for all the studied points, we 
find the highest sensitivity to 1) 5^ b R , 5f b L and 5^; 2) 5 R b R the next, 3) 5^ t R the 
next to next; and 4) the lowest sensitivity is to 5 RL and 5 RR . Consequently, these 
two later will be the less constrained ones. 

- Comparing the predictions with the experimental data: 

In this case, the experimental allowed 3<T ex p band is very narrow, and all the 
central predictions at 5* Y = 0, i.e. for MFV scenarios, lay indeed outside this 
band. However, if we assume again that our predictions suffer of a similar large 
uncertainty as the SM prediction, given in (I45p . then the MFV predictions are all 
compatible with data except for SPS4. It is also obvious from this figure that the 
predictions within NMFV, as compared to the MFV ones, lie quite generically 
outside the experimental allowed band, except for the above commented deltas 
with low sensitivity. 

- Intervals of 5^ Y allowed by data: 

We consider again, that a given bf Y value is allowed by AM Bs data if the predicted 
interval AM Ba ± A theo (AMg s ), intersects the experimental band. It corresponds 
to adding linearly the experimental uncertainty and the theoretical uncertainty. 
Given the present controversy on the realistic size of the theoretical error in the 
estimates of A theo (AMsJ in the MSSM (see, for instance, [73]), we choose a very 
conservative value for the theoretical error in our estimates, considerably larger 
than the SM value in (gBJ), of A thco (AM Bs ) = 51 x 10~ 10 MeV. This implies that 
a predicted mass difference in the interval 

63 x 10~ 10 < AM Bs (MeV) < 168.6 x 1(T 10 , (48) 

is regarded as allowed. 

The allowed intervals for the deltas that are obtained with this requirement are 
collected in table HJ As we have already commented, the restrictions on the b- 
sector parameters from AM Ba are very strong, and in consequence, there are 
narrow intervals allowed for, 5^ R , S R b L , and 5 23 ■ In the case of 5 R b R there are 
indeed sequences of very narrow allowed intervals, which in some cases reduce to 
just a single allowed value. The parameters that show a largest allowed interval, 
with sizeable |<5^ y |, larger than 0(0.1), are 5 RR , S RL and 5^ R . 
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Figure 2: Sensitivity to the NMFV deltas in BR(B -> X s j) for the SPSX points of table HJ 
The experimental allowed 3a area is the horizontal colored band. The SM prediction and 
the theory uncertainty A thco (BR(5 — > Xg^y)) (red bar) is displayed on the right axis. 
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Figure 3: Sensitivity to the NMFV deltas in BR(B S -> for the SPSX points of table HJ 

The experimental allowed region by the 95% CL bound is the horizontal colored area. The 
SM prediction and the theory uncertainty A theo (BR(i? s — > fi + fi~)) (red bar) is displayed on 
the right axis. 
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experimental allowed 3cr exp area is the horizontal colored band. The SM prediction and the 
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Total Allowed 5j) Y Intervals 



We finally summarize in table |3] the total allowed intervals for all the NMFV deltas, S^ Y , 
where now we have required compatibility with the present data of the three considered B 
observables, BK(B — > Xg'y), BK(B S — > and AM# s . It is obvious, from the previous 

discussion, that the most restrictive observables are BR(I? — > X s 7)and AM^ s , leading to a 
pattern of allowed delta intervals which is clearly the intersect of their two corresponding 
intervals. The main conclusion from this table is that, except for SPS4 (the point SPS4 is 
practically excluded), the NMFV deltas in the top-sector can be sizeable \5^. Y \ larger than 
O(0. 1) and still compatible with B data. In particular, 5 RL , and 8 RR are the less constrained 
parameters, and to a lesser extent also S^ t R . The parameters on the bottom-sector are, in 
contrast, quite constrained. The most tightly constrained are clearly 5^ h R and S R b L , specially 
the first one with just some singular allowed values: either positive and of the order of 
3 — 5 x 10 -2 , or negative and with a small size of the order of —7 x 1CT 3 ; for the second the 
limits are around 2 x 10~ 2 for both positive and negative values. 5 R b R is the less constrained 
parameter in the bottom sector, with larger allowed intervals of \8 RR \ ^ 0.4 — 0.9 depending 
on the scenario. 

All SPS points are defined with a positive \i value. We have checked the effect of switching 
the sign of \x. While the numerical results are changing somewhat, no qualitative change can 
be observed. Consequently, confining ourselves to positive /i does not constitute a general 
restriction of our analysis. Similar observations are made in the Higgs-sector analysis below. 

The intervals allowed by B data that we have presented above will be of interest for the 
following study in this work, where we will next explore the size of the radiative corrections 
to the MSSM Higgs masses within these NMFV-MSSM scenarios and we will require com- 
patibility with B data. In the final analysis of these corrections, we will use the constraints 
from B data as extracted from two non-vanishing deltas. As expected, these constraints 
vary significantly respect to the ones with just one non-vanishing delta. 
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BR(B -> X al ) 


BR{B S -4 H+/J.-) 


AM fl , 




SPSla 
SPSlb 
SPS2 
SPS3 
SPS4 
SPSS 


(-0.51:-0.43) (-0.034:0.083) 
(-0.33:-0.27) (-0.014:0.062) 

(-0.43:0.34) (0.90:0.92) 
(-0.73:-0.65) (-0.083:0.12) 
(-0.14:-0.11) (0.0069:0.034) 
(-0.26:0.50) 


(-0.53:0.92) 
(-0.014:0.16) 

(-0.99:0.99) 

(-0.90:0.97) 
(0.028:0.055) 

(-0.60:0.57) 


(-0.73:-0.65) (-0.41:0.55) (0.73:0.79) 
(-0.090:-0.069) (-0.021:0.097) (0.14:0.17) 
(-0.37:0.37) 
(-0.86:-0.79) (-0.56:0.66) (0.83:0.89) 
(-0.0069) (0.021:0.055) (0.076) 
(-0.37:0.39) 




SPSla 
SPSlb 


(-0.89:-0.86) (-0.12:-0.097) (-0.062:0.28) 
(-0.083:0.36) 


(-0.89:0.89) 
(-0.44:0.67) 


(-0.89:0.89) 
(-0.67:0.67) 


SLR 
°ct 


SPS2 
SPS3 


(-0.46:0.46) 
(-0.43:0.61) 


(-0.46:0.46) 
(-0.68:0.68) 


(-0.46:0.46) 
(-0.68:0.68) 




SPS4 


(-0.61:-0.51) (0.041:0.23) 


excluded 


(-0.39:-0.021) (0.74:0.77) 




SPS5 


(-0.27:0.58) 


(-0.59:0.61) 


(-0.59:0.61) 




SPSla 
SPSlb 
SPS2 
SPS3 
SPS4 
SPS5 


(0) (0.034) 
(-0.0069:0) (0.048:0.055) 
(-0.0069:0) (0.048:0.055) 
(-0.0069:0) (0.048:0.055) 
(-0.0069) (0.034) 
(-0.0069:0) (0.041) 


(-0.60:0.60) 
(-0.43:0.54) 
(-0.48:0.48) 
(-0.61:0.61) 

(0.49) 
(-0.71:0.71) 


(-0.076:0.076) 

(-0.15:0.14) 

(-0.19:0.19) 

(-0.12:0.12) 
(-0.29:-0.24) (-0.10:-0.014) (0.12:0.18) 
(-0.090:0.090) 




SPSla 


(-0.84:0.84) 


(-0.84:0.84) 


(-0.84:0.84) 




SPSlb 


(-0.63:0.63) 


(-0.63:0.63) 


(-0.63:0.63) 




SPS2 
SPS3 


(-0.39:0.39) 
(-0.63:0.63) 


(-0.39:0.39) 
(-0.63:0.63) 


(-0.39:0.39) 
(-0.63:0.63) 




SPS4 


excluded 


excluded 


(-0.72:-0.21) (0.21:0.72) 




SPS5 


(-0.53:0.53) 


(-0.53:0.53) 


(-0.53:0.53) 


t>3, 1 


SPSla 
SPSlb 
SPS2 
SPS3 


(-0.014:0.014) 
(-0.021:0.021) 
(-0.014:0.014) 
(-0.021:0.021) 


(-0.71:0.71) 
(-0.58:0.58) 
(-0.55:0.55) 
(-0.63:0.63) 


(-0.069:0.069) 
(-0.14:0.14) 
(-0.17:0.17) 
(-0.11:0.11) 




SPS4 


(-0.021:-0.014)(0. 014:0.021) 


excluded 


(-0.21:-0.17) (0.16:0.21) 




SPSS 


(-0.014:0.014) 


(-0.72:0.72) 


(-0.083:0.083) 




SPSla 


(-0.93:-0.67) (-0.64:0.93) 


(-0.93:0.93) 


(-0.93:0.93) 




SPSlb 


(-0.93:-0.61) (-0.56:0.90) 


(-0.95:0.94) 


(-0.98:0.98) 




SPS2 
SPS3 
SPS4 


(-1.0:0.99) 
(-0.97:0.97) 
excluded 


(-1.0:0.99) 
(-0.97:0.97) 
excluded 


(-1.0:0.99) 
(-0.98:0.97) 
(-0.85:-0.22) (0.22:0.85) 




SPS5 


(-0.60:0.60) 


(-0.60:0.60) 


(-0.60:0.60) 




SPSla 


(-0.65:0.68) 


(-0.96:0.96) 


(-0.91.--0.90) (-0.86:-0.80) (-0.41:0.41) (0.81:0.86) (0.90:0.91) 




SPSlb 
SPS2 
SPS3 


(-0.71:0.74) 
(-0.99:0.99) 
(-0.98:0.98) 


(-0.73:0.98) 
(-0.99:0.99) 
(-0.98:0.98) 


(-0.94:-0.92) (-0.83:0.88) (0.93:0.94) 

(-0.99) (-0.39:0.39) (0.99) 
(-0.94:-0.93) (-0.88:0.88) (0.93:0.94) 




SPS4 
SPSS 


(-0.45:-0.18) (0.19:0.46) 
(-0.77:0.80) 


excluded 
(-0.97:0.97) 


(-0.80:-0.028) (0.461:0.71) (0.86:0.91) (0.94:0.95) 
(-0.92) (-0.87:-0.78) (-0.51:0.51) (0.78:0.87) (0.92) 







Total allowed intervals 




SPSla 


(-0.034:0.083) 




SPSlb 


(-0.014:0.062) 


S.LL 
"23 


SPS2 
SPS3 


(-0.37:0.34) 
(-0.083:0.12) 




SPS4 


(0.028:0.034) 




SPS5 


(-0.26:0.39) 




SPSla 


(-0.89:-0.86) (-0.12:-0.097) (-0.062:0.28) 




SPSlb 


(-0.083:0.36) 




SPS2 


(-0.46:0.46) 


SPS3 


(-0.43:0.61) 




SPS4 


excluded 




SPS5 


(-0.27:0.58) 




SPSla 
SPSlb 


(0)(0.034) 
(-0.0069:0) (0.048:0.055) 


S.LR 
°sb 


SPS2 
SPS3 


(-0.0069:0) (0.048:0.055) 
(-0.0069:0) (0.048:0.055) 




SPS4 


excluded 




SPS5 


(-0.0069:0) (0.041) 




SPSla 


(-0.84:0.84) 




SPSlb 


(-0.63:0.63) 


xRL 

°ct 


SPS2 
SPS3 


(-0.39:0.39) 
(-0.63:0.63) 




SPS4 


excluded 




SPS5 


(-0.53:0.53) 




SPSla 


(-0.014:0.014) 




SPSlb 


(-0.021:0.021) 


SjRL 
°sb 


SPS2 
SPS3 


(-0.014:0.014) 
(-0.021:0.021) 




SPS4 


excluded 




SPS5 


(-0.014:0.014) 




SPSla 


(-0.93:-0.67) (-0.64:0.93) 




SPSlb 


(-0.93:-0.61) (-0.56:0.90) 


xRR 
°ct 


SPS2 
SPS3 


(-1.0:0.99) 
(-0.97:0.97) 




SPS4 


excluded 




SPS5 


(-0.60:0.60) 




SPSla 


(-0.41:0.41) 




SPSlb 


(-0.71:0.74) 


zRR 
°sb 


SPS2 
SPS3 


(-0.99) (-0.39:0.39) (0.99) 
(-0.94:-0.93) (-0.88:0.88) (0.93:0.94) 




SPS4 


excluded 




SPS5 


(-0.51:0.51) (0.78:0.80) 



Table 3: Total allowed delta intervals by BR(B X s j), BR(B S fi + > _ ) and AM Ba . 
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4 Radiative corrections to MSSM Higgs masses within 
NMFV scenarios 



In this section we present our computation of the one-loop radiative corrections to MSSM 
Higgs boson masses within the NMFV scenarios. For completeness and definiteness, we 
first shortly review the relevant features of the MSSM Higgs sector at tree-level. Then we 
summarize the main one-loop renormalization issues that are involved in the computation 
and finally we present the analytical results for the one-loop corrected Higgs masses. 



4.1 The Higgs boson sector at tree- level 

Contrary to the SM, in the MSSM two Higgs doublets are required. The Higgs potential [74 



V 



m{\Utf + mi\H 2 \' - m{ 2 {e ab UlU b 2 + h.c.) 

+\{gl + gl) [W -\u 2 \ 2 ] 2 + \g 2 2 \'h\u 2 \ 



(49) 



contains mi,m 2 , as soft SUSY breaking parameters; g 2 , g\ are the SU{2) and U (1) gauge 
couplings, and e\ 2 = —1. 

The doublet fields H\ and H 2 are decomposed in the following way: 



Hi 



- <xS) 



«2 + ^2 + *X2J 



(50) 



The potential (l4"9"j) can be described with the help of two independent parameters (besides 
g 1 and g 2 ): tan/3 = v 2 /vx and M\ = — m 2 2 (tan/3 + cot/3), where Ma is the mass of the 
CP-odd Higgs boson A. 

The diagonalization of the bilinear part of the Higgs potential, i.e. of the Higgs mass 
matrices, is performed via the orthogonal transformations 



H 

h 
G 
A 

G ± 



cos a sin a 

- sin a cos a 
cos (3 sin j3 

- sin j3 cos (3 
cos (3 sin j3 

- sin {3 cos (3 



X? 



The mixing angle a is determined through 

-(Ml + M|)sin/3cos/3 



a = arctan 



Mf cos 2 P + M\ sin 2 (3-ml Mee 
One gets the following Higgs spectrum: 

2 neutral bosons, CV = +1 : 



7T 

- - < a < 



(51) 
(52) 
(53) 

(54) 



h,H 
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1 neutral boson, CP = — 1 
2 charged bosons 
3 unphysical Goldstone bosons 



.4 

H + ,H~ 



(55) 



At tree level the mass matrix of the neutral CP-even Higgs bosons is given in the <f>\- 
basis in terms of M z , Ma, and tan/3 by 



M 



2, tree 
Higgs 



m 



m 



•1 <!>■> 
2 



<f>i<fa ""<p2 

A\ sin 2 P + M% cos 2 - (Ml + Mf ) sin (3 cos (3 
(M\ + M 2 )sin/3cos/3 M| cos 2 /3 + M| sin 2 (3 



(56) 



which by diagonalization according to Eq. ( 15T|) yields the tree-level Higgs boson masses 



M 2,tree 
1V1 Higgs 



///.tree 







ft, tree 



(57) 



where 



/J tree 



Ml + M l ± 



(Mi 



The charged Higgs boson mass is given by 



M 2 ) 2 -4M 2 Micos 2 20 



(58 



(59) 



The masses of the gauge bosons are given in analogy to the SM: 



' 1l 2 i 2\ 



Ml 



M~ = 0. 



(60) 



4.2 The Higgs boson sector at one-loop 

In order to calculate one-loop corrections to the Higgs boson masses, the renormalized Higgs 
boson self-energies are needed. Here we follow the procedure used in Refs. [HUES] (and 
references therein) and review it for completeness. The parameters appearing in the Higgs 
potential, see Eq. (1491) . are renormalized as follows: 



Mf Mf + 8M 2 Z , 
M^^M^ + 5M^, 

Mffiggs — Mffiggs + ^^Higgs; 



T h T h + 5T h , 
T H ^T H + 5T H , 
tan/3 — >• tan/3 (1 + <5tan/3 ). 



(61) 



Mg iggs denotes the tree- level Higgs boson mass matrix given in Eq. (156]) . and T# are the 
tree-level tadpoles, i.e. the terms linear in h and H in the Higgs potential. 

The field renormalization matrices of both Higgs multiplets can be set up symmetrically, 




\SZhh 



Uz h 



tin 



1 + \5Z HH 




(62) 



24 



For the mass counter term matrices we use the definitions 

5m 2 h Sm hH 



6M£, m = I " 1 . (63) 

Hlsgs 1 5m 2 hH 6m 2 H l 

The renormalized self-energies, X(p 2 ), can now be expressed through the unrenormalized 
self-energies, £(p 2 ), the field renormalization constants and the mass counter terms. This 
reads for the CP-even part, 

t hh (p 2 ) = Z hh (p 2 ) + 5Z hh (p 2 - m 2 h>tiee ) - 6m 2 h , (64a) 
Z hH (p 2 ) = Z hH (p 2 ) + 5Z hH {p 2 - |(m 2 trec + m| itre J) - 5m 2 hH , (64b) 
T lHH (p 2 ) = Y, HH (p 2 ) + 5Z HH (p 2 - m 2 HtTCC ) - 5m 2 H . (64c) 

Inserting the renormalization transformation into the Higgs mass terms leads to ex- 
pressions for their counter terms which consequently depend on the other counter terms 
introduced in ( 16"TT) . 

For the CP-even part of the Higgs sectors, these counter terms are: 

5m 2 h = 5M 2 A cos 2 (a - 0) + 5M 2 Z sin 2 (a + /3) (65a) 
+ 2m/s w c w (^Th cos(a - 0) sin 2 (a - 0) + 5T h sin(a - 0)(1 + cos 2 (a - /?))) 
+ 5tan/3 sin cos f3 (M 2 A sin 2 (a - /3) + M§ sin 2(a + 0)), 
5m 2 hH = \(5M 2 A sm2(a- (3) - 5Mf sin2(a + /3)) (65b) 

+ mk^( 5T * sin3 (« - P) ~ 6T * cos3 ( a " ft) 

- <5tan/3 sin/3 cos/3 (M A cos 2 (a - /3) + Mf cos2(a + /3)), 

= 5M| sin 2 (a - /3) + 5Mf cos 2 (a + /3) (65c) 
~ 2Mzs wCv ( ST h cos ( a ~ P)Q- + sin 2 (a - 0)) + <5T h sin(a - /3) cos 2 (a - /3)) 

— 5tan/3 sin /3 cos /3 (M A sin 2(a — /3) + M| sin 2(a + /3)) . 

For the field renormalization we chose to give each Higgs doublet one renormalization 
constant, 

Hi -+ (1 + ISZnJHu U 2 ^{1 + \5Z H2 )U 2 . (66) 

This leads to the following expressions for the various field renormalization constants in 
Eq. (J62D: 



SZ hh = sm 2 a SZ ni + cos 2 a SZ H2 , (67a) 
5ZhH = sin a cos a (5Zu 2 — dZ^), (67b) 
SZ HH = cos 2 a SZ Hl + sin 2 a 5Z H2 . (67c) 

The counter term for tan /3 can be expressed in terms of the vacuum expectation values as 

aan/3 = \ (6Z n2 - SZ Hl ) + - , (68) 
2 v 2 v\ 
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where the 8vi are the renormalization constants of the vf. 

vt ->■ (1 + SZ-Hj) (vi + Svi) , v 2 -> (1 + SZ H2 ) (v 2 + Sv 2 ) 



(69) 



Similarly for the charged Higgs sector, the renormalized self-energy is expressed in terms of 
the unrenormalized one and the corresponding counter-terms as: 



where, 
and, 



£h-h+ (p 2 ) = ^h-h+ (p 2 ) + SZ H - H + (p 2 - m 2 H± trcc ) - 5m 2 H± , 
Sm 2 H ± = 8M 2 A + 8M^ 



(70) 

(71) 
(72) 



5Z H - H + = sin 2 (3 5Z Hl + cos 2 (3 5Z n2 . 

The renormalization conditions are fixed by an appropriate renormalization scheme. For 
the mass counter terms on-shell conditions are used, resulting in: 



5M 2 = Re£ zz (Mf), bM 2 , = ReE w (M 2 ,), SM 2 A = Re£ AA (M 2 ). 



(73) 



For the gauge bosons £ denotes the transverse part of the self-energy. Since the tadpole 
coefficients are chosen to vanish in all orders, their counter terms follow from TUrn. + 
$T{ h ,H} = 0: 



ST, = —Th 



5T 



-T 



H 



(74) 



For the remaining renormalization constants for <5tan/3, 5Z%^ and 5Zu 2 the most convenient 
choice is a DR renormalization of <5tan/3, 5Z^ and 5Zu 2 , 



DR. 



5Z Hl = 5Z Hi 
6Z H2 = 5Z™ 



1 div 



<5tan/3 = --{8Z %2 - 5Z ni 



[Re T! HH | a=0 ] 
[Re T! hh | a=0 ] , 

5ttm(5 m 



(75a) 
(75b) 

(75c) 



The corresponding renormalization scale, /i^, is set to /ij^ = m t in all numerical evalua- 
tions. 

Finally, in the Feynman diagrammatic (FD) approach that we are following here, the 
higher-order corrected CP-even Higgs boson masses are derived by finding the poles of the 
(h, if )-propagator matrix. The inverse of this matrix is given by 



(Am f 



P 



m 



tree + ^HH(P 2 ] 



t hH ( P 2 ) 



Z hH {p 2 ) p 2 -ml tvcc + E hh (p 2 



(76) 



Determining the poles of the matrix Amggs in Eq. ( 1761) is equivalent to solving the equation 

1 2 



P ~ m h,tree + ^hh(p 



P 



m 



H±rcc 



(P 2 



^hH(p 



0. 



(77) 



Similarly, in the case of the charged Higgs sector, the corrected Higgs mass is derived by the 
position of the pole in the charged Higgs propagator, which is defined by: 



2 2 
P ~ m H±,trcc 



t H - H+ {p 2 ) = 0. 



(7? 
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4.3 Analytical results of Higgs mass corrections in NMFV-SUSY 



Following the previously detailed prescription for the computation of the one-loop corrected 
Higgs boson masses, one finds the analytical results for these masses in terms of the renor- 
malized self-energies which, in turn, are written in terms of the unrenormalized self-energies 
and tadpoles. To shorten the presentation of these analytical results, it is convenient to 
report just on these unrenormalized self-energies and tadpoles. 

The relevant one-loop corrections have been evaluated with the help of FeynArts [38] 
and FormCalc [39]. For completeness the new Feynman rules included in the model file are 
listed in the Appendix A. All the results for the unrenormalized self-energies and tadpoles 
are collected in Appendix B. We have shown explicitly just the relevant contributions for the 
present study of the radiative corrections to the Higgs boson masses within NMFV scenarios, 
namely, the one-loop contributions from quarks and squarks. The corresponding generic 
Feynman-diagrams for the Higgs bosons self-energies, gauge boson self-energy diagrams and 
tadpole diagrams are collected in the Fig. [15] in Appendix B. It should also be noticed that 
the contributions from the squarks are the only ones that differ from the usual ones of the 
MSSM with MFV. It should be noted also that the corrections from flavor mixing, which 
are the subject of our interest here, are implicit in both the Vckm, and in the values of the 
rotation matrices, R u , R d , and the squark masses, m^, m^. (i = 1, .., 6) that appear in these 
formulas of the unrenormalized self-energies and tadpoles and that have been introduced in 
section [21 

Finally, it is worth mentioning that we have checked the finiteness in our analytical results 
for the renormalized Higgs self-energies. It is obviously expected, but it is not a trivial check 
in the present scenarios with three generations of quarks and squarks and with flavor mixing. 
We have also checked that the analytical results of the self-energies in Appendix B agree with 
the formulas in FeynHiggs [T6l[T7T[2Tl[22] . Each one of the terms contained in the Appendix 
B was compared with the corresponding term in FeynHiggs. During this process and the 
check of the finiteness, discrepancies were found with the charged Higgs part of FeynHiggs, 
leading to an updated version of the co 



5 Numerical analysis 

In this section we present our numerical results for the radiative corrections to the Higgs 
boson masses from from flavor mixing within NMFV-SUSY scenarios. Since all one-loop 
corrections in the present NMFV scenario are common to the MSSM except for the correc- 
tions from squarks, which depend on the bf Y values, we will focus just on the results of 
these corrections as a function of the flavor mixing parameters, and present the differences 
with respect to the predictions within the MSSM. Correspondingly, we define: 

Am,(5f) = < MPV (C) - ™f SM > <f> = h,H,H±, (79) 

where m^ MFV (5$ Y ) and m^ ssu have been calculated at the one- loop level. It should be noted 
that m 1 } MFV (5f j Y = 0) = m ^ SSM and, therefore, by construction, Am^^f = 0) = 0, and 
Am^ gives the size of the one-loop NMFV contributions to m^. The numerical calculation of 



7 We especially thank T. Halm for his efforts put into this update. 
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m NMFV^xy-j anc j m MSSM |- )een d one with (the updated version of) FeynHiggs p!6 | [T7 | [21 ] 
[22J , which solves the eqs. [77J and [75] for finding the positions of the poles of the propagator 
matrix. Previous results for ArrihiS^) can be found in [20] . 

5.1 Am0 versus one d* Y ^ 

We show in figs. [5] [6J and [7J our numerical results for Arrih, ArriH and Am#±, respectively, 
as functions of the seven considered flavor changing deltas, 5^, S^ R , 5^ R , S RL , 5 R b L , S RR and 
5 R b R , which we vary in the interval — 1 < bf Y < 1. In these plots we have chosen the same 
six SPS points of tabled], as for the previous study of constraints from B physics in El We do 
not take the experimental bounds into account here, since we just want to show the general 
behavior of the masses with the deltas. The experimental bounds will be taken into account 
in the next subsection. As before we have checked the impact of switching the sign of /i and 
found a small quantitative but no qualitative effect. 

The main conclusions from these figures are the following: 

- General features: 

All mass corrections, Arrih, Amu and Am#±, are symmetric bf Y — > —5* Y , as ex- 
pected. This feature is obviously different than in the previous plots of the B observ- 
ables. The sign of the mass corrections can be both positive and negative, depending 
on the particular delta value. The size of the Higgs mass corrections, can be very large 
in some 5^ Y ^ regions, reaching values even larger than 10 GeV in some cases, at 
the central region with not very large delta values, \5^ Y \ < 0.5. In fact, the restrictions 
from B physics in this central region is crucial to get a reliable estimate of these effects. 

For low tan (3, where the restrictions from B physics to the deltas are less severe, 
the Higgs mass corrections are specially relevant. Particularly, Am/, turns out to be 
negative and large for tan/3 = 5 (SPS5) for all deltas, except 8 R b R . For instance, at 
\S^ Y \ ~ 0.5, the mass correction Arrih for SPS5 is negative and > 5 GeV in all flavor 
changing deltas except 5 R b R where the correction is negligible. In the case of Amu and 
Am#± the size of the correction at low tan j3 is smaller, < 2 GeV in the central region, 
except for 5^ R and 5 R b L that can also generate large corrections > 5 GeV. 

In the cases with large tan/3 (SPS4 and SPSlb), we also find large mass corrections 
but, as already said, they are much more limited by B constraints. In particular, for 
SPS4 all deltas are excluded, except for a very narrow window in <5^ 3 L (see tabled]). 

In the cases with moderate tan/3 = 10 (SPSla, SPS2 and SPS3), we find large correc- 
tions | Arrih | ~ 5 GeV in the central region of S^ b R , S R b L , S^ t R and S RL . The other Higgs 
bosons get large corrections |Am#|, |Am#| 5 GeV in the deltas central region only 
for 8g b R and 5f b L . 

- Sensitivity to the various deltas: 

We find very strong sensitivity in the three mass corrections Arrih, Amu and Am#±, 
to 8g R and S R b L for all the seven considered SPS points. 

In the case of Arrih there is also an important sensitivity to 5^ R and S RL in all the 
considered points. The strong sensitivity to LR and RL parameters can be understood 
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due to the relevance of the A-terms in these Higgs mass corrections. It can be noticed 
in the Feynman rules (i.e. see the coupling of two squarks and one/two Higgs bosons 
in Appendix A) that the A-terms enter directly into the couplings, and in some cases, 
as in the couplings of down-type squarks to the CP-odd Higgs boson, enhanced by 
tan /3. Therefore, considering the relationship between the A-terms and these LR and 
RL parameters as is shown in eq. |2U the strong sensitivity to these parameters can 
be understood. A similar strong sensitivity to 5^ R in Am^ has been found in [26] . 

In SPS5 there is a noticeable sensitivity to all deltas except 5^ R . In other points, the 
effects of , S RR on Am/, are only appreciated at the large delta region, close to ±1. 
For instance, in SPS2, Am^ = —5 GeV for S RR = ±1. 

In the case of Am#, apart from S^ R and 6f b L , there is only noticeable sensitivity to 
other deltas in SPS5. The same comment applies to Am#±. 
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Figure 5: Sensitivity to the NMFV deltas in Am/, for the SPSX points of tabled) 
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5.2 Am/, versus two d^ Y ^ 



13 



Our previous results on the Higgs mass corrections show that the corrections to the lightest 
Higgs mass Am/, are negative in many of the studied cases and can be very large for some 
regions of the flavor changing deltas which are still allowed by present B data. These negative 
and large mass corrections, can lead to a prediction for the corrected one-loop mass in these 



kind of NMFV-SUSY scenarios, 



m NMFV ^ 



m 



MSSM 
h 



+ Am h , which are indeed too low and 



already excluded by present data [6j[7]. Therefore, interestingly, the study of these mass 
corrections can be conclusive in the setting of additional restrictions on the size of some 
flavor changing deltas which otherwise are not bounded from present B data. 

In order to explore further the size of these 'dangerous' mass corrections, we have com- 
puted numerically the size of Am/, as a function of two non- vanishing deltas and have looked 
for areas in these two dimensional plots that are allowed by B data. We show in Figs. [9], [10l 
[TTj [T2l fl3| and HH the numerical results of the Ami, contour-lines in two dimensional plots, 
(5^,6*- Y ), for the respective points BFP, SPS2, SPS3, SPS5, VHeavyS and HeavySLightH 
of table dJ 



BR(B -> X sl ) 



BR(B -> X s7 ) 





Figure 8: Legend for plots of Higgs mass corrections varying two deltas simultaneously 
displayed in figs. |9l [TQl [TT1 [121 El] and HH Each colored area represents the disallowed 
region by the specified observable/s inside each box. A white area placed at the central 
regions of the mentioned figures represents a region allowed by the three B observables. A 
white area placed outside the colored areas represent regions of the parameter space that 
generate negative squared masses. These problematic points are consequently not shown in 
our plots, as we did in the previous plots. The discontinuous lines in those figures represent 
the contour lines for the B observables corresponding to the maximum and minimum allowed 
values: dash-dot-dash for the upper bound of BR(£> — > X s 7)(eq. (|46|) ). dot-dash-dot for the 
lower bound of BR(_B — > X s 7)(eq. (H6l)). dashed line for the upper bound of AMg s (eq. 
( H8"|) ). a sequential three dotted line for the lower bound of AM Bg (eq. (H8"|) ). and a dotted 
line for the upper bound of BR(B S — > /i + /i~)(eq. pT]) ). 
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Figure 9: Am/, (GeV) contour lines from our two deltas analysis for BFP. The color code 
for the allowed/disallowed by B data is given in figJHJ 
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Figure 10: Am/, (GeV) contour lines from our two deltas analysis for SPS3. The color code 
for the allowed/disallowed by B data is given in figJHJ 
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Figure 11: Am^ (GeV) contour lines from our two deltas analysis for SPS2. The color code 
for the allowed/disallowed by B data is given in figJHJ 
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Figure 12: Am^ (GeV) contour lines from our two deltas analysis for SPS5. The color code 
for the allowed/disallowed by B data is given in figJHJ 
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Figure 14: Am/, (GeV) contour lines from our two deltas analysis for HeavySLightH. The 
color code for the allowed/disallowed areas by B data is given in figJHJ 
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We have chosen in all plots as one of these non- vanishing deltas mainly because of two 
reasons. First, because it is one of the most frequently studied flavor changing parameters 
in the literature and, therefore, a convenient reference parameter. Second, because there are 
several well motivated SUSY scenarios, where this parameter gets the largest value, as we 
explained in section P2J 

In these two-dimensional figures we have included the allowed/disallowed by B data 
areas that have been found by following the procedure explained in section |3l and the 
allowed intervals are given in eqs. (l46j) . ( I47j) and (HE])- The color code explaining the meaning 
of each colored area and the codes for the discontinuous lines are given in figM Contour 
lines corresponding to mass corrections above 60 GeV or below -60 GeV have not been 
represented. In several scenarios the plots involving 5^ R ' R show a seemingly abrupt behavior 
for \8 S R,RL \ ^ 0.3, corresponding to extremely large (one-loop) corrections to m^. In general, 
in the case of very large one-loop corrections, in order to get a more stable result further 
higher order corrections would be required, as it is known from the higher-order corrections 
to m h in the MFV case (see, e.g., Ref. [H]). However, we cannot explore this possibility 
here. On the other hand, in order to understand the behavior of as a function of 5^ R ' RL 
a simple analytical formula would have to be extracted from the general result. However, 
this is beyond the scope of our paper. 

The main conclusions from these two dimensional figures are summarized in the following: 

The points that have been chosen in these plots are quite representative of all the different 
patterns found. The plots for SPSla (not shown here) manifest similar patterns as those of 
SPS3. The plots for SPSlb (not shown here) manifest similar patterns as those of BFP. The 
plots for SPS4 are not included because they do not manifest any allowed areas by B data. 

The largest mass corrections Am/, found, being allowed by B data occur in plots (f^^, S^ t R ) 
and (623 , 5 RL ). This applies to all the studied points. They can be as large as (—50, —100) 
GeV at 8^ or 5 RL close to the upper and lower horizontal bands in these plots where 5^ R 
or S RL are close to ±0.5. Again these large corrections from the LR and RL parameters are 
due to the A-terms, as we explained at the end of section [5711 Comparing the different plots, 
it can be seen that the size of the allowed area by the B data (the white area inside of the 
colored regions) can be easily understood basically in terms of tan/3, and the heaviness of 
the SUSY and Higgs spectra. Generically, the plots with largest allowed regions and with 
largest Higgs mass corrections correspond to scenarios with low tan /3 = 5 and heavy spectra. 
Consequently, the cases of VHeavyS and HeavySLightH are the most interesting ones, ex- 
hibiting very large radiative corrections, resulting from the heavy SUSY spectra. In the case 
of HeavySLightH the large corrections are not only found for Am/,, but also, though to a 
lesser extent, for the other Higgs bosons, Am# and Am#± (not shown here). Consequently, 
in this scenario the deltas will be very restricted by the mass bounds, especially by m^. 

There are also important corrections in the allowed areas of the two dimensional plots 
of (623 , 8 RR ) for some points, particularly for SPS5 (and to a lesser extent for SPS2). Here 
the corrections can be as large as -50 GeV in the upper and lower parts, i.e. for 5 RR close 
to ±0.5. In the case of SPS2 they can be up to -2 GeV for this same region. 

As for the remaining two-dimensional plots they do not show relevant allowed areas where 
the mass corrections are interestingly large. 
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6 Conclusions 



In this paper we have analyzed the one-loop corrections to the Higgs boson masses in the 
MSSM with Non-Minimal Flavor Violation. We assume the flavor violation is being gener- 
ated from the hypothesis of general flavor mixing in the squark mass matrices, and these are 
parametrized by a complete set of S^ Y (X, Y = L,R; i,j —t,c,u or b, s, d). 

In the first step of the analysis we scanned over the NMFV parameters, contrasting 
them with the experimental bounds on BR(I? — > X s j), BK(B S — > /i + /i~) and AM# S . We 
take into account the most up-to-date evaluations in the NMFV MSSM for BR(5 — > X s ~/), 
BK(B S — > and AMg s , as included in the BPHYSICS subroutine of the SuFla code [3T] . 

For the evaluation of AMb s we have added the one-loop gluino boxes [32] which are 
known to be very relevant in the context of NMFV scenarios [331435] . We have estimated 
the size of these corrections and compared them with the other relevant contributions from 
chargino boxes and double Higgs penguins for all values of tan (5 for the first time. And we 
have concluded that gluino boxes dominate for moderate and low tan (3 < 20 which is the 
interesting range for the present work. In the final part of the B physics analysis, we have 
evaluated in one-dimensional scans which intervals for the 5^ Y are still allowed in certain 
benchmark scenarios based on the SPS points. 

In the second step we analyzed the one-loop contributions of NMFV to the MSSM Higgs 
boson masses, focusing on the parameter space still allowed by the experimental flavor con- 
straints and by current limits from Higgs boson searches. Here two relevant 6f Y were varied 
simultaneously, thus enlarging the allowed range for these parameters. We found large cor- 
rections, mainly for the low tan case, up to several tens of GeV for rrih and somewhat 
smaller corrections for m# and m#±. These corrections are specially relevant in the case of 
the light MSSM Higgs boson since they can be negative and up to two orders of magnitude 
larger than the anticipated LHC precision. Consequently, these corrections must be taken 
into account in any Higgs boson analysis in the NMFV MSSM framework. Conversely, in 
the case of a Higgs boson mass measurement these corrections might be used to set further 
constraints on bf^ ' . The present work clearly indicates that the flavor mixing parameters 
5c t R and 5 RL are severely constrained by the present bounds on the lightest Higgs boson mass 
within the NMFV-MSSM scenarios. 
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Appendix A 

We list the new Feynman rules of the NMFV scenario that are involved in the present 
computation. The corresponding couplings to the Higgs boson H are obtained from the 
ones listened here for the lightest Higgs boson h by replacing 

C a — > S a ; S a -> — C a ; S a+ p — > —C a+ p ; C2a ~^ ~C2a (80) 

The notation used for the formulas is the following: s x = sinx; c x = cosx; s w = 
sm9w] c w = cosOw = j^rf', tp = tan/3. 

1. Couplings of two squarks and one/two Higgs bosons 



Ui 

h° 

< 



h° 



H+ 



■ -< 



-< 



di 



-< 



■Em=i 6Mff I Sw , W {Qc a c w ml k - M w M z s a+ ,s,(3 - 4s 2 J) 

+3c w R u j3+l {A u kl c a m Uk + S k i m Uk fi*s a )} 

+3c w Rf l (Af* k c a m Ul + S k i ni Uk fis a ) } } 



! n ElUi , Mw Zs^ {*h (6wr& - M w M z s a+p c p {3 - 2sD) 

dj N " +3c w Rj j3+l (Af l s a m dk + S k i m dk ffc a ) j 

+Rf,l +k {hiRis+i (Qs a c w m 2 dk - 2M w M z s a+f) cpsl) 
+3c w R d jl (A d * k s a m dl + S k i m dk fic a ) J | 



J2k,l=l y/2MwSwtR {J2n=l An*k m «»^CKM + m u k /^CKM^) 

<*i ^ +m dl m Uk V^ l M Ri 3+ i(l + 1/)} 

+ R i,k {"^3+^/3 (J2n=l A n,l m 4^CKM^ + m d l /^*Vckm) 

+^ckm Rj,i R fc " tf,(M$y8 2p - m 2 



^1,1=1 y/2Mws*tf> \ R iUl { R lk l P (En=l A tl m dn^CKM^ + m d t A'^CKm) 

+m dl m Uk V^ M Rl 3+k (l + t p 2 )} 
+V£kuRh (m 2 Uk - t p {M^s 2p - m 2 t p ))}} 
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-J2k,l=l 2Mws^t p { R i,3+k R ll ( A l,k m ui + S k im Uk fitp) 
-Yfk,l=l 2Mw Sw { R i,3+k R j,l ( A f,k m dltp + Sium dh fi) 

- R U R U + i (Ki m ^ H + m d k m*) } 

+2Rl 3+k Rj* +k (3c 2 a clml k - 2c 2a M^s}sl) } 
~ ELl 12M^444 {^fc^fc ( 6s a C w m 4 + c 2a M^c^(3 - 24 )) 

+2 R t3+k R ll + k (3^c>l + c 2q m^44) } 
"X)fc=i i2mJ54^ { R i,k R j,k ( 3c l m t k - M w s U 3 ~ 4s «)) 

+ R f, 3+k R U +k (3c>^ - 2M£^4) } 

- ELi vmiks^ { R l* R U ( 6c w< - ^M^(3 - 44)) 

+2 R h +k R Ts + k (3c>^ fc - 2c w M^ty w )} 
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A di 

- el H^fe ( 6c >i^ + c ^( 3 - 2s w)) 

/ \ \/y W W ^ 

A ° dj +2Ris + k R ls+k ( 3c >l*l + 0^4) } 

if" Ui 



,, /3 

2 M i o^2 



if- 4 



X _ V^J ief J od* pd (aST^ 6 ™2 T/nfe T/*rii „2 

^ ' " 2^k,l=l 12A/^ S 2 C 2^ \ n i,k n j,l 1° 2^n=l m «„ CKM ''CKM^ 

H+ dj +c 2 „5 kl M^tl(l-4cl)) 

+26 kl R% +k Rj >3+l t} (3clt}m 2 dk + c w M^s\ 



)} 



3. Couplings of two squarks and one/two gauge bosons 



w~ 



di 



dj 



u 



dj 



ELi idk (^RtURh+A - m R U^ ~ 4s w)) (p + p ) 



ELi *fc (2^3 +fc ^ 3+fc 4 - i#i&(3 - 24)) (p + p ) 



E^=i^^c^^(p+p), 
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di y 



U; 



$2jfe=l 18c 2 s 2 (-^ffc-^i,fc(3 2s^) 2 + 4:Rf^ + kRj, 3+k S t ) 



!/,; 



E3 te 2 p« * pi 
fe=l 2sl n 'i,k rL j,k 



di 



E3 ie 2 pd* pd „ 
fc=l 2s 2 £ H,k rL j,k 



d J 
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Appendix B 

All the following Feynman diagrams have been calculated using FeynArts 3.5 [38] and Form- 
Calc 6.0 [39]. The notation used here is the same as in Appendix A. Furthermore we use 
the functions [76] 

* a r 21 f fJ. 4 ~ D d D k 1 

A \m 2 \ = / n— — (81) 

16tt L J J (2tt) d k - m 2 K J 

" B tf, ml ml] . / l^^—^-^—^ (82) 

p> Bl [pi m\, m i] = f ^-^- 2 ifr ^— 5 ^ (83) 

J (27r) [k 2 - mf J [(A; + p) - m l 2 \ 

The generic diagrams have been ordered according to its topologies, and the particles 
involved in the internal loops (quarks q or squarks q). The diagrams can be found in figlTBl 
The complete self-energy can be expressed as a sum of three parts: 

V V 2| J _L_ V 2 *? _L V 1 ? V V 2< ? _L_ V 2| 3 _L_ V 1 ? T -1 T"3 _L T"? (QA\ 

~ + ^u' + - z-vv + ^VV + ^VV 1 <t>- 1 4> + 1 4> 

where 0, 0' = /i, H, A, H and V = W, Z. All the self-energies £ correspond to X (p 2 ). 
The self-energies for H are obtained by the replacements of eqlEO] on the results of h: 
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Figure 15: Different topologies for E<^/, Hyv, 
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The self-energies for Th are obtained using the replacements of eq. (150|l on the results 
of T h . 
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